Irreducible finite-size effects in surface free energies from crystal-nucleation data. 
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In this Letter we report a simulation study in which we compare the solid-liquid interfacial free 
energy of NaCl at coexistence, with the value that follows from the height of the homogeneous 
nucleation barrier. We find that the two estimates differ by more than 100%. Similar, although 
smaller discrepancies are found for crystals of hard-sphere colloids and of Lennard- Jones ("argon") 
particles. We consider a variety of possible causes for this discrepancy and are forced to conclude 
that it is due to a finite-size effect that cannot be corrected for by any simple thermodynamic 
procedure. Importantly, we find that the surface free energies that follow from real nucleation 
experiments should be subject to large finite size effects. Taking this in to account, we obtain 
quantitative agreement between the simulation data and the surface free energy of NaCl that follows 
from nucleation experiments. Our finding suggests that most published solid-liquid surface free 
energies derived from nucleation experiments will have to be revised. 
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The study of homogeneous crystal nucleation is of in- 
terest because it provides information about the pathway 
by which crystalline order emerges from the disordered 
parent phase. However, such experiments are also of con- 
siderable practical importance, as they are used to esti- 
mate the magnitude of the solid-liquid interfacial free en- 
ergy. Classical Nucleation Theory (CNT, see e.g. [l[) pro- 
vides the route by which experimental nucleation rates 
are related to surface free energies. CNT relates the 
number of crystal nuclei that form per second per cu- 
bic meter (denoted by R) to AG cri t, the height of the 
free-energy barrier that has to be crossed to nucleate a 
crystaki? = k e - AG cru/k B T ^ jj ere K j s a kinetic prefac- 
tor, Tis the absolute temperature and ks is Boltzmann's 
constant. CNT predicts the following expression for the 

height of the nucleation barrier: AG cr u — c -^r^ji > 
where 7ls is the liquid-solid surface free energy per unit 
area, Ap is the difference in chemical potential between 
the solid and the supercooled liquid, and ps is the num- 
ber density of the crystalline phase, c is a constant that 
depends on the shape of the cluster, e.g. c = 167r/3 
for a spherical crystal nucleus. As the nucleation rate 
depends exponentially on AG cr it, the rate is a very sen- 
sitive function of the surface free-energy density jls- A 
crucial assumption underlying CNT is that the bulk and 
surface properties of a small crystal nucleus are the same 
as those of a macroscopic crystal. However, it has been 
long realized that this assumption is questionable, as a 
critical crystal nucleus often contains only a few hundred 
molecules. Indeed, in his review on crystal nucleation, 
Kelton writes: "...while the precise meaning of [tls] is 
uncertain, it constitutes a parameter that can be deter- 
mined for each element and profitably used to make pre- 



dictions of the nucleation behavior". In other words: the 
7ls determined from nucleation experiments can only be 
used to predict the outcome of other nucleation experi- 
ments, thus severely limiting the predictive value of CNT. 
More microscopic theories such as density-functional the- 
ory (DFT) the Cahn-Hilliard approach (CH) 0, 
or the phase-field formalism (PF) [5| can, and have been 
used to improve upon CNT. Yet, the question remains 
whether the widely-used CNT can be reformulated in 
such a way that it correctly describes the properties of 
small clusters yet, at te same time, reproduces the cor- 
rect, macroscopic surface free energy. 

Increasingly accurate simulation techniques allow us to 
probe both the free energy of small nuclei and the surface 
free energies of planar crystal-liquid interfaces. A case in 
point is the system NaCl in contact with its melt. Ref. @ 
reported the surface free energy of a NaCl [100] interface 
in contact with the coexisting liquid phase: 7ls ~ 36 ± 6 
mJ m -2 . The effective surface free energy that follows 
from the NaCl crystal-nucleation barrier at 800 K was 
reported in ref. [7]: 7ls = 80 ± 1 mJ m~ 2 (assuming 
that the nucleus has a cubic shape). In addition, nucle- 
ation experiments at 905K 0] provide an experimental 
estimate of 7ls ~ 68 mJ m~ 2 . 

Another example of a large difference between 7ls de- 
rived from the nucleation barrier and from coexistence 
data comes from hard-sphere colloids: a comparison of 
simulations at coexistence Q and in the supersaturated 
regime [l(| indicate that the value of 7ls estimated on 
the basis of the nucleation barrier is some 30% larger than 
the value for a planar interface at coexistence. A simi- 
lar discrepancy (O(20%)) exists between the surface free 
energy for the planar interface and the crystal nucleus 
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of the (truncated and (force-)shifted) Lennard- Jones po- 
tential [IlT - [l3j ]. It is clearly of considerable interest to 
understand the origin of the discrepancy between the 
nucleation data and the results for 7ls at coexistence, 
as this might facilitate the interpretation and analysis of 
experimental nucleation data. 

In this Letter, we report a systematic study of the 
finite-size effects in the surface free energy of NaCl crys- 
tals in contact with their melt. We chose this system 
because it shows the largest discrepancy of all exam- 
ples listed above. As in refs. @, 0, EH we use the Tosi- 
Fumi rigid-ion-pair interaction potential 15| to model the 
inter-ionic interactions in NaCl. We note at the outset 
that it is imprecise to speak of the surface free-energy 
density of a small crystallite, as the value of 7ls depends 
on the choice of the dividing surface (equimolar dividing 
surface, equi-enthalpy dividing surface, surface of tension 
etc. - see ref. [16]). For flat interfaces, the correspond- 
ing surface free energies are all the same, but this is not 
the case for strongly curved surfaces. The surface free 
energy that enters into CNT is the one associated with 
the surface of tension 17| . One property of the surface of 
tension is that it is, to lowest order, independent of the 
choice of the dividing surface. We use this property to 
determine 7ls associated with the surface of tension. To 
facilitate the comparison with the data of ref. [6] that re- 
fer to a flat interface at coexistence, we deduce 7ls from 
the size-dependence of the free energy of a small crystal- 
lite at coexistence. At coexistence, there is no difference 
in chemical potential between the liquid and the bulk 
solid, hence the excess free energy of a small crystallite 
is entirely due to its surface. 

All simulations were carried out at the coexistence tem- 
perature T M =(1060±10) K. T M reported in [I1EI1- The 
melting temperature of Tosi-Fumi NaCl is close to the 
experimental value: T^ p —1074 K As a first step, 

we determine the dependence of the free energy of small 
NaCl crystallites on the number of ions in the crystal. 
For this part of the calculation, we make use of umbrella 
sampling [20| at constant N, P and T. These simulations 
yield the excess free energy of the largest crystalline clus- 
ter in the system as a function of the number of parti- 
cles in that cluster. We use a geometrical criterion (see 
ref. [7j) to distinguish crystalline from liquid- like parti- 
cles. We can then deduce the surface free-energy density 

using 



7ls(A0 



AG(N) 
C 



p 2 J 3 N~ 2 / 3 . 



(1) 



where N denotes the number of crystalline particles and 
C is a geometrical constant that has a value 6 for a cubic 
nucleus and (3671") 1 / 3 for a sphere. Although there are 
strong fluctuations in the shape of a small NaCl crystal- 
lite in contact with its melt, its average shape is fairly 
cubic (see figure [T]) . Of course, we need not assume a 
priori that the cluster is cubic: we can use the average 




FIG. 1: Average shape of a NaCl crystallite consisting of 140 
solid-like particles. As the shape of such a small cluster fluc- 
tuates, we obtain the average shape by superimposing a large 
number of instantaneous configurations of the same mass, fix- 
ing their center of mass, and the orientation of the crystal 
axes. We average over all 48 symmetry-related orientations. 
The surface is defined as the set of points where the average 
density equals the average of the solid and liquid densities. 
Only the crystalline particles inside this surface are shown. 



cluster shape from figure [T] to perform a Wulff construc- 
tion (see e.g. [2l[) that yields the variation of the surface 
free energy with orientation. Assuming that the surface 
free energy of the [100] equals the macroscopic value, we 
can the compute the average 7ls of the cluster. Leav- 
ing apart the question whether a Wulff construction is 
at all meaningful for clusters containing O(10 2 ) parti- 
cles, we note that this procedure yields (7ls) ~ 40 mJ 
mT 2 , which is within 10% of the value expected for a 
perfect cube. In what follows, we will therefore assume 
that small NaCl crystals have the same cubic morphology 
as macroscopic crystals. From Figure [2] we cannot yet 




FIG. 2: Free energy barriers (/3AG) as a function of N g at 
coexistence Tm- The error bars in /3AG are comparable to 
the size of the symbols. 

deduce the surface free energy because there is no a pri- 
ori reason to assume that the surface of this geometrical 
cluster has any thermodynamic meaning. We know, how- 
ever, that in the thermodynamic limit, the ratio of N g to 
N, the "thermodynamic" number of atoms in the crys- 
tal, should approach 1. We therefore make the ansatz: 



N = (Ng /3 +a) 3 , where a is an adjustable parameter that 
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remains to be determined. To find the number of atoms 
within the surface of tension, we choose a value of a that 
minimizes the variation of 7ls with the size of the clus- 
ter. This analysis leads to a value of a pa 0. Figure ^ 
shows that, over range of cluster sizes studied, the result- 
ing value of 7ls is indeed almost independent of N for 
all but the smallest clusters. More interestingly, we find 
that the resulting value of 7ls is very close to the value 
7ls ~ 80m J m~ 2 that follows from the analysis of the 
nucleation barrier at 800K [3] (see figured]). Moreover, a 
similar analysis at 800 K, leads to the same estimate of 
7ls- The internal consistency between the values of 7ls 




FIG. 3: Interfacial free energy as a function of N computed 
at the surface of tension, assuming a cubic cluster with a=0 

derived from the nucleation barrier and from the surface 
of tension would be encouraging, were it not for the fact 
that it does nothing to resolve the discrepancy with the 
value of 36 mJ m -2 found for a flat interface. Choosing 
another conventional dividing surface (e.g. the equimo- 
lar or the equi-enthalpy surface) [23| only makes matters 
worse: in both cases we find a negative value of a that 
results in an even larger value of 7ls that is, moreover, 
strongly cluster-size dependent. Hence, we conclude that 
the discrepancy between the properties of a small crystal 
and a flat interface cannot be corrected for by choosing 
a better definition of the location of the solid-liquid in- 
terface. 

The above discussion suggests that the conventional 
version of CNT cannot account for the observed discrep- 
ancy between the surface free energy of a flat interface 
and that of a small crystallite. However, even within a 
thermodynamic approach, one can introduce corrections 
to classical nucleation theory that would change the ap- 
parent value of the surface free energy. One such correc- 
tion takes into account that the crystal nucleus is com- 
pressible and that the surface free energy depends on the 
density of the crystal. To estimate the magnitude of this 
effect, we extend the analysis of Mullins [17( to obtain: 

(7Ls(p s )/7Ls(p s (0)) 3 



where r is the ratio between the barrier height in the 
case of compressible nuclei (with density p s ), compared 
to that for incompressible clusters (with density p s (0))- 
B denotes the bulk modulus of the crystal and e the 
elastic strain, compared to that of a solid at the same 
chemical potential as that of the parent liquid. From our 
simulations, we find that the density at the center of the 
crystal nucleus is some 6% lower than the reference value. 
Even with this rather extreme estimate of the strain in 
the nucleus, we find a compressibility-correction to the 
apparent value of 7ls that is no more than 10%. Hence, 
we conclude that compressibility effects cannot account 
for the observed discrepancy. 

Thus far, we have not considered the effect of edges 
and vertices on the surface free energy of a small cluster. 
This effect is certainly non- negligible. If, for instance, 
we consider a cubic NaCl crystal in vacuum at T = OK, 
both the line energy of the edges and the vertex energy 
of the corners can be determined directly. The energy of 
an NaCl cube can be written as: 



e B £ + 6e s r + l2e E e + 8e c 



(2) 



(1 - A Ws /(2B) - \Be^/(p s ^))' 



where e is the total internal energy per particle, e E , the 
energy per particle in a bulk crystal, es the energy of 
a particle belonging to the surface, e E the energy of a 
particle belonging to an edge, ec the energy of a particle 
belonging to a corner of the cube, and £ is the number 
of atoms per edge. Computing this energy for a crystal 
of 64, 216 and 512 atoms [Hj], we find that e E /e s = 0.22 
and ec/e s = 1.2. The effect of these edge and vertex 
contributions is to increase the apparent surface energy 
by 13% for a crystal of 216 particles. Of course, these 
numbers do not apply to a hot NaCl crystal in contact 
with its melt and it is not at even obvious how to define 
the various terms in that not only the magnitude 

but even the sign of e E and ec depend on the precise 
choice of the dividing surface. This means that, within 
the macroscopic framework imposed by CNT, we cannot 
reliably estimate the edge and corner contributions to the 
surface free energy. 

We are therefore forced to conclude that the large ap- 
parent value of 7ls of small crystallites is due to a finite- 
size effect that is not easily accounted for by within a 
"thermodynamic" theory. Rather, the free energy of small 
clusters must be computed using a molecular approach, 
either theoretically (as in DFT [2, 3],CH 4], or PF [|) 
or numerically, as illustrated in the present work. In the 
present paper, and in ref. Q, we computed the free energy 
of relatively small clusters (up to iV=200). However, un- 
der the experimental conditions for crystal nucleation of 
NaCl (T=905K), the critical nucleus is expected to con- 
tain 0(6 x 10 2 ) particles. Calculations for larger clusters 
would be feasible, but expensive. We therefore make the 
We therefore make the Tolman "ansatz" that the lead- 
ing correction to surface free energy is proportional to 
l/R e , where R c is the radius of the critical nucleus. As 
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1/Rc ~ A/i, we assume that the variation in 7ls is of 
the form: 7ls(A^)= 7£g ox + bA/i. We can determine 
b from the simulation data of refs [fl 0]- Inserting the 
value Afj, = 0.3kT at T=905K, we predict that under 
the condition of the nucleation experiments of ref. Q, 
the effective value of 7ls should be 67 mJ m~ 2 , in al- 
most embarrassing agreement with the experimental data 
(7ls=68 mJ m -2 ). Although this good agreement is al- 
most certainly fortuitous, it does support our conjecture 
that the surface free energies measured in nucleation ex- 
periments are subject to very large finite size corrections 
(in this case: more than 80%). If we take this strong 
A/z-dependence of 7ls seriously, it would mean that for 
strongly faceted crystals (although not for NaCl), the nu- 
cleation barrier could start to rise again at large super- 
saturations. This should be experimentally observable, 
as it would lead to an increase in the final crystallite size 
in fully crystallized samples j24[. Interestingly - but we 
do not know if it is really relevant - the final crystallite 
size in hard-sphere crystallization suddenly grows as the 
concentration is increased beyond a volume fraction of 
58%. If the barrier is a monotonically decreasing func- 
tion of the volume fraction, this should not happen. 

In summary, our study of the free energy of NaCl crys- 
tallites indicates that the surface free energy is subject 
to large finite size corrections that cannot be accounted 
for within a thermodynamic theory. Based on the small 
number of examples where the relevant simulation data 
are available (NaCl, Lennard Jones, hard spheres), we 
speculate that the finite size effects are most pronounced 
for strongly faceted crystals, such as NaCl. The present 
results support the suggestion by Kelton that the large 
number of published surface free energies that are based 
on nucleation data are of little use to predict macroscopic 
surface free energies. We stress that, in addition to nucle- 
ation studies, there are other, more reliable, experimental 
routes to determine solid-liquid surface free energies. An 
example is the grain-boundary groove method 25] . How- 
ever, such experiments are challenging, especially for ma- 
terials that have anisotropic surface free energies. Our 
work highlights the fact that, if nucleation studies are 
used to estimate solid-liquid surface free energies, the 
analysis cannot be based on CNT but must make us of 
one of the more accurate, microscopic theories for crystal 
nucleation that properly account for the fact that crystal 
nuclei are far from macroscopic. 
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